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Executive  Summary 

Classification  and  Discrimination  of  Sources 
with  Time-Varying  Frequency  and  Spatial  Spectra 

ONR  Grant  No.  N000 14-98-1-01 76 
Moeness  Amin  (PI) 

This  report  presents  the  results  of  the  research  performed  under  ONR  grant  number  N00014- 
98-1-0176  over  the  period  of  October  1st,  2004  to  September  30th,  2005.  The  research  team 
working  on  this  project  consists  of  Prof.  Moeness  Amin  (PI),  Prof.  Yimin  Zhang  (Research 
Professor),  Mr.  Baha  Obeidat  and  Mr.  Pawan  Setlur  (Graduate  Students).  We  have  also 
collaborated  with  Dr.  Thayananthan  Thayaparan  (Defense  Research  and  Development,  Canada), 
Prof.  Ahmad  Hoorfar  (Villanova  University),  and  Prof.  Kehu  Yang  (Xidian  University,  China). 

The  research  objective  in  FY05  evolved  around  the  moving  target  tracking  using  time- 
frequency  signal  representations.  In  addition  to  achieving  this  objective,  the  research  has 
progressed  on  three  other  fronts,  namely.  Sensor  configuration  in  polarized  antenna  arrays, 
DOA  and  polarization  estimation  for  wideband  sources,  and  Micro-Doppler  experiments.  One 
book  chapter,  two  journal  articles,  and  seven  conference  papers  have  been  generated  from  our 
research  efforts  in  FY05.  Below  is  the  summary  of  the  research  accomplishments  in  each  of  the 
individual  areas.  Copies  of  the  principle  publications  are  included  in  the  report. 

1.  Spatial  Polarimetric  Time-Frequency  Distributions  for  Moving  Target  Tracking 

We  have  introduced  the  spatial  polarimetric  time-frequency  distributions  (SPTFDs)  as  a 
platform  for  processing  polarized  nonstationary  signals  incident  on  multiple  dual-polarized 
double-feed  antennas.  Based  on  this  platform,  we  have  developed  the  polarimetric  time- 
frequency  MUSIC  (PTF-MUSIC)  methods  for  direction-of-arrival  (DOA)  estimation  of 
nonstationary  sources  with  distinct  polarization  characteristics,  and  have  examined  the  feasibility 
of  the  PTF-MUSIC  methods  for  tracking  moving  sources  with  time-varying  polarization 
characteristics.  We  have  demonstrated  the  significance  of  polarization  diversity  in  challenging 
direction  finding  problems,  where  the  sources  are  closely  spaced,  and  discuss  important  issues 
relevant  to  utilization  of  polarization  diversity.  The  SPTFD  has  the  capability  of  incorporating 
both  the  instantaneous  polarization  information  and  time-frequency  signatures  of  the  different 
sources  in  the  field  of  view.  The  incorporation  of  specific  time-frequency  points  or  regions, 
where  one  or  more  signals  reside,  enhances  signal-to-noise  ratio  (SNR)  and  allows  source 
discrimination  and  source  elimination.  This,  in  turn,  leads  to  DOA  performance  improvement 
and  reductions  in  the  required  number  of  array  sensors.  The  PTF-MUSIC  significantly 
outperforms  the  existing  time-frequency  MUSIC,  polarimetric  MUSIC  and  conventional  MUSIC 
direction  finding  techniques. 
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2.  DOA  and  Polarization  Estimation  for  Wideband  Sources 

We  have  presented  an  approach  to  direction-of-arrival  (DOA)  and  polarization  estimations  of 
far-field  wideband  signals  using  an  array  of  cross-polarized  dual-feed  sensors.  The  wideband 
signals  are  decomposed  into  a  set  of  narrowband  signals.  Coherent  signal  subspace  processing  is 
then  applied  to  dimension-extended  source  signal  correlation  matrices  at  different  narrowband 
frequencies,  leading  to  the  construction  of  single  signal  subspace.  The  utilization  of  polarization 
diversity  through  the  use  of  a  dual-feed  cross-polarized  array  to  process  the  polarized  signal 
arrivals  lends  an  additional  degree  of  freedom,  leading  to  performance  improvement  over  the 
case  where  a  single-polarization  array  is  used. 

3.  Sensor  Configurations  in  Nonstationary  Array  Processing 

.  We  have  developed  the  spatio-polarimetric  correlation  as  a  framework  to  evaluate  the 
direction-of-arrival  estimation  performance  of  different  receiver  configurations  employing 
double-feed  dual-polarized  array  as  well  as  arrays  of  single-feed  single-polarized  elements. 
Single-  and  dual-polarized  array  configurations  are  often  used  in  imaging  and  radar  applications. 
The  spatio-polarimetric  correlation  coefficients  are  derived  for  several  array  reconfigurations  and 
used  to  provide  a  configuration-based  performance  comparison  based  on  spatial  resolution  and 
grating  lobes. 

4.  Micro-Doppler  Experiments 

We  have  conducted  several  important  indoor  experiments  for  detecting  micro-Doppler 
signatures  of  some  commonly  encountered  indoor  targets.  Targets  in  real  life  have  vibrating  and 
rotating  parts  in  addition  to  the  translations  endured  by  the  target.  The  target  itself  might  either 
be  rotating  or  vibrating  with  no  translation,  a  motion  is  described  mechanically  as  simple- 
harmonic  motion  (SHM).  Such  nonlinear  motions  then  induce  additional  frequency  modulations 
around  the  Doppler,  this  is  termed  micro-Doppler  effect.  Micro-DOppler  has  many  applications 
for  target  recognition  and  classification.  Experiments  were  designed  to  image  targets  which 
induce  micro-Doppler  like  signatures.  The  targets  were  primarily  a  rotating  fan  and  a  toy  train 
moving  on  a  circular  track.  The  experiments  were  conducted  using  CW  (Continuous  Wave) 
radars.  These  radars  measure  only  the  real  component  of  target  returns.  Time-frequency 
techniques  were  used  to  depict  the  instantaneous  frequency  of  the  target.  The  objective  of  the 
experiments  was  to  compute  the  frequency  of  vibrations  or  rotations.  In  this  report,  we  introduce 
the  problem  statement,  and  propose  a  model  for  Micro-Doppler.  Experimental  setup  is  described, 
and  simulation  and  experimental  results  are  also  presented. 
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Polarimetric  Time-Frequency  MUSIC  for  Direction  Finding 
of  Moving  Sources  with  Time-Varying  Polarizations 

Yimin  Zhang,  Baha  A.  Obeidat,  Moeness  G.  Amin 

Abstract:  The  spatial  polarimetric  time-frequency  distributions  (SPTFDs)  have  been 
introduced  as  a  platform  for  processing  polarized  nonstationary  signals  incident  on  multiple 
dual-polarized  double-feed  antennas.  Based  on  this  platform,  the  polarimetric  time-frequency 
MUSIC  (PTF-MUSIC)  methods  have  been  developed  for  direction-of-arrival  (DOA) 
estimation  of  nonstationary  sources  with  distinct  polarization  characteristics.  In  this  paper, 
we  examine  the  feasibility  of  the  PTF-MUSIC  methods  for  tracking  of  moving  sources  with 
time- varying  polarization  characteristics.  We  demonstrate  the  significance  of  polarization 
diversity  in  challenging  direction  finding  problems,  where  the  sources  are  closely  spaced,  and 
discuss  important  issues  relevant  to  utilization  of  polarization  diversity.  The  SPTFD  has  the 
capability  of  incorporating  both  the  instantaneous  polarization  information  and  time- 
frequency  signatures  of  the  different  sources  in  the  field  of  view.  The  incorporation  of 
specific  time-frequency  points  or  regions,  where  one  or  more  signals  reside,  enhances  signal- 
to-noise  ratio  (SNR)  and  allows  source  discrimination  and  source  elimination.  This,  in  turn, 
leads  to  DOA  performance  improvement  and  reductions  in  the  required  number  of  array 
sensors.  The  PTF-MUSIC  significantly  outperforms  the  existing  time-frequency  MUSIC, 
polarimetric  MUSIC  and  conventional  MUSIC  direction  finding  techniques. 
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1.  Introduction 


Time-frequency  distributions  (TFDs)  have  been  used  for  nonstationary  signal  analysis 
and  synthesis  in  various  areas,  including  speech,  biomedicine,  automotive  industry,  and 
machine  monitoring  [1,  2,  3].  Over  the  past  few  years,  the  spatial  dimension  has  been 
incorporated,  along  with  the  time  and  frequency  variables,  into  quadratic  and  higher-order 
TFDs,  and  led  to  the  introduction  of  spatial  time-frequency  distributions  (STFDs)  as  a 
powerful  tool  for  nonstationary  array  signal  processing  [4,  5,  6].  The  relationship  between  the 
TFDs  of  the  sensor  data  to  the  TFDs  of  the  individual  source  waveforms  is  defined  by  the 
steering,  or  the  mixing,  matrix,  and  was  found  to  be  similar  to  that  encountered  in  the 
traditional  data  covariance  matrix  approach  to  array  processing.  This  similarity  has  allowed 
subspace-based  estimation  methods  to  be  applied  to  time-frequency  signal  representations, 
utilizing  the  source  instantaneous  frequency  for  direction  finding.  It  has  been  shown  that  the 
MUSIC  [7,  8],  maximum  likelihood  [9],  and  ESPRIT  [10]  techniques  based  on  SPTFDs 
outperform  their  counterparts  that  are  based  on  covariance  matrices,  when  the  sources  are  of 
nonstationary  temporal  characteristics. 

Polarization  and  polarization  diversities,  on  the  other  hand,  are  commonly  used  in 
wireless  communications  and  various  types  of  radar  systems  [11,  12].  Antenna  and  target 
polarization  properties  are  widely  employed  in  remote  sensing  and  synthetic  apertureradar 
(SAR)  applications  [13,  14,  15].  Airborne  and  spacebome  platforms  as  well  as 
meteorological  radars  include  polarization  information  [16,  17].  Additionally,  polarization 
plays  an  effective  role  in  identification  of  targets  in  the  presence  of  clutter  [18,  19],  and  has 
also  been  incorporated  in  antenna  arrays  to  improve  signal  parameter  estimations,  including 
DO  A  and  time-of-arrival  (TOA)  [20, 21]. 
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To  combine  time-frequency  (t-f)  signal  representations  and  polarimetric  signal 
processing,  we  have  introduced  the  spatial  polarimetric  time-frequency  distributions 
(SPTFDs)  for  double-feed  dual-polarized  arrays,  which  allow  the  utilization  of  both  the  t-f 
and  polarization  signatures  of  the  sources  in  the  field  of  view  [22,  23].  We  then  used  the 
SPTFD  to  develop  the  polarimetric  time-frequency  MUSIC  (PTF-MUSIC)  technique,  which  . 
is  formulated  based  on  the  source  t-f  and  polarization  properties.  This  new  technique  is 
employed  for  DOA  estimation  of  polarized  nonstationary  signals  and  shown  to  outperform 
the  MUSIC  techniques  that  only  incorporate  either  the  t-f  or  the  polarimetric  source 
characteristics  [23].  The  application  of  the  SPTFDs  to  an  ESPRIT-like  method  was 
introduced  in  [22,24], 

In  this  paper,  we  examine  the  feasibility  of  the  PTF-MUSIC  methods  for  tracking  of 
moving  sources  with  time-varying  polarization  characteristics.  We  demonstrate  the 
significance  of  polarization  diversity  in  direction  finding  problems  in  which  the  sources  are 
closely  spaced,  and  discuss  the  necessary  treatment  of  the  data  samples  to  effectively  utilize 
the  polarization  diversity.  The  SPTFD  incorporates  the  instantaneous  polarization 
information  at  selected  auto-term  t-f  points  where  the  energy  of  desired  source  signals  is 
localized.  The  selection  of  t-f  regions  in  constructing  the  SPTFD  matrices  not  only  reduces 
noise,  but  also  enables  source  exclusion  from  post-eigendecomposition  and  subspace 
estimation.  Reduction  of  the  number  of  sources  reduces  the  required  number  of  array  sensors 
and  lead  to  performance  enhancement.  As  a  result,  PTF-MUSIC  outperforms  existing 
polarimetric  MUSIC  technique. 

This  paper  is  organized  as  follows.  Section  0  discusses  the  signal  model  and  briefly 

reviews  TFDs  and  STFDs.  Section  0  considers  dual-polarized  antenna  arrays  and  introduces 
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the  concept  of  spatial  polarimetric  time-frequency  distributions  (SPTFDs).  The  PTF-MUSIC 
algorithm  is  proposed  in  Section  0.  Section  0  considers  the  spatio-polarimetric  correlation 
which  is  a  useful  parameter  to  measure  the  closeness  between  two  sources  in  the  joint  spatial 
and  polarimetric  domains.  Section  0  considers  DOA  tracking  of  moving  sources  with  time- 
varying  polarization  characteristics.  Computer  simulations,  demonstrating  the  effectiveness 
of  the  proposed  methods,  are  provided  in  Section  0. 

Throughout  this  paper,  vectors  and  matrices  are  denoted  by  lower-case  bold  and  upper¬ 
case  letters,  respectively.  The  operation  £[•]  denotes  expectation,  (•)*  denotes  complex 
conjugate,  (-)T  denotes  transpose,  and  (-)H  denotes  conjugate  transpose  (Hermitian).  We 
use  to  denote  polarization,  ||  •  ||  vector  norm,  0  denotes  the  Kronecker  product  operator, 
and  o  denotes  the  Hadamard  product  operator. 

2.  Signal  Model 

2.1  Time-Frequency  Distributions  (TFDs) 

The  Cohen's  class  of  TFDs  of  a  signal  x(t)  is  defined  as  [1] 

DxxitJ)  =  jJW*  ~  u,t)x(u  +  ^)x*  (u  -  ~)e"j24rdudT ,  (1) 

where  t  and / represent  the  time  and  frequency  indexes,  respectively,  0(t,r)  is  the  t-f  kernel, 
and  t  is  the  time-lag  variable.  In  this  paper,  all  the  integrals  are  from  -  oo  to  oo . 

The  cross-term  TFD  of  two  signals  x,  (t)  and  xk  (t)  is  defined  by 

DXlxk  /)  =  -  u,  t)X'  (u  +  l)x;  (u  -  ^)e~jWTdudr .  (2) 
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2.2  Spatial  Time-Frequency  Distributions  (STFDs) 

The  STFDs  have  been  developed  for  single-polarized  antenna  arrays  [4,  5,  6]. 
Consider  n  narrowband  signals  impinging  on  an  array  of  m  single-polarized  antenna 
elements.  The  following  linear  data  model  is  assumed, 

x(0  =  y(0  +  “(0  =  A(O)s(0  +  n(/),  (3) 

where  the  mxn  matrix  A(<D)  =  [a(^,),a(^2 a(^n)]  is  the  mixing  matrix  that  holds  the 
spatial  information  of  the  n  signals,  0  =  [^,,^2,...,^n]  and  a(^,.)  is  the  spatial  signature  for 
source  /.  Each  element  of  the  «xl  vector  s(t)  =  [s,  (/),  s2  (t),  ...,sn  (f )] 7  is  a  mono-component 
signal.  Due  to  the  signal  mixing  occurring  at  each  sensor,  the  elements  of  the  mxl  sensor 
data  vector  x(/)  are  multi-component  signals.  n(t)  is  an  mxl  additive  noise  vector  that 
consists  of  independent  zero-mean,  white  and  Gaussian  distributed  processes. 

The  STFD  of  a  data  vector  x(t)  is  expressed  as  [4] 

DIX(t,/)  =  jfy(t-u,r)x(u  +  ~)xH (u-~)e~j2,lfTdudT ,  (4) 

where  the  (/, k)  th  element  of  DX1  (t,  f)  is  given  by  Eq.  (2)  for  i, k  =  1,2, The  noise-free 
STFD  is  obtained  by  substituting  Eq.  (3)  in  Eq.  (4),  resulting  in 

Dxx(r,/)  =  A(0)Dss(r,/)AK(0),  (5) 

where  Dss(t,/)  is  the  TFD  matrix  of  s(t)  which  consists  of  auto-  and  cross-source  TFDs. 
With  the  presence  of  the  noise,  which  is  uncorrelated  with  the  signals,  the  expected  value  of 
Dm(T/)  yields 
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£[D„(',/)]  =  AmE[J)ss(t,f)]AH(0)  +  a2l , 


(6) 


where  a1  is  the  noise  power  and  I  is  the  identity  matrix. 

Equation  (6)  relates  the  STFD  matrix  to  the  source  TFD  matrix.  It  is  similar  to  the 
formula  that  is  commonly  used  in  narrowband  array  processing  problems,  relating  the  source 
covariance  matrix  to  the  sensor  spatial  covariance  matrix.  It  is  clear,  therefore,  that  the  two 
subspaces  spanned  by  the  principle  eigenvectors  of  Dxx  (t,f)  and  the  columns  of  A(<D)  are 
identical.  By  constructing  the  STFD  matrix  from  the  t-f  points  of  highly  localized  signal 
energy,  the  corresponding  signal  and  noise  subspace  estimates  become  more  robust  to  noise 
than  their  counterparts  obtained  using  the  data  covariance  matrix  [5,  8,  9].  Further,  source 
elimination,  rendered  through  the  selection  of  specific  t-f  regions,  improves  DOA 
estimations. 

3.  Spatial  Polarimetric  Time-Frequency  Distributions  (SPTFDs) 

3.1  Polarimetric  Modeling 

For  a  far-field  transverse  electromagnetic  (TEM)  wave  incident  on  the  array,  shown  in 
Figure  1,  the  electric  field  can  be  described  as 

E(t)  =  Ee(t)0  +  E ,{tyf 

=  [Eg  (t)  COS(0 )  cos(^)  -  Ej  (0  sin(^)]x  (7) 

+  [Eg  (0  cos(<9)  sin(^)  -  E*  ( t )  cos(^)]j>  +  Ee  (t)  sin  (0)z 

where  ^  and  6  are,  respectively,  the  spherical  unit  vectors  along  the  azimuth  and  elevation 
angles  ^  and  6 ,  viewed  from  the  source.  The  unit  vectors  x ,  y  and  z  are  defined  along  the 
x ,  y  and  z  directions,  respectively.  For  simplicity  and  without  loss  of  generality,  it  is 
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assumed  that  the  source  signal  is  in  the  x—y  plane,  whereas  the  array  is  located  in  the  y—z 

A 

.  plane.  Accordingly,  6  =  90° ,  0  =  -z,  and 

E(t)  =  -Et(t)sm(</>)x  +  Ej(t)cos(0)y  +  E()(t)z.  (8) 

We  denote  s(t)  as  the  source  signal  magnitude  measured  at  the  receiver  reference  sensor, 
with  polarization  angle  y  e  [0,90°  ],  and  polarization  phase  difference  rj  e  (-180°  ,180°] .  The 
source’s  vertical  and  horizontal  polarization  components,  s[v](t)  and  s[h](t),  can  then  be 
expressed  in  terms  of  the  respective  spherical  fields,  Ee(t)  and  EM),  as 

Ee(t)  =  s[v](t)  =  s(t)cos(y),  Ej(t)  =  s[h](t)  =  s(t)sm(y)eJn .  (9) 

A  signal  is  referred  to  as  linearly  polarized  if  r]  =  0°  or  rj  =  180°.  Substituting  Eq.  (9)  into 
Eq.  (8)  results  in 

E(/)  =  5(t)[-cos(y)sin(^)x  +  cos(^)sin(y)e7,7j)  +  cos(y)z].  ( 10) 

Now  we  consider  that  n  signals  impinge  on  the  array,  consisting  of  m  dual-polarized 
antennas.  The  vertical  and  horizontal  components  of  the  i  th  source  are  expressed  as 

5/Vl(0  =  si  (0cos(/i)  =  cnS'  (t),  sjA](t)  =  Sj  (t)sin(yi)eJr,i  =caSj(t),  (  11) 

where  the  parameters  cn  =cos(yi)  and  cj2  =sin (yi)eJ’’1  denote  the  vertical  and  horizontal 
polarization  coefficients.  The  signal  received  at  the  /  th  dual-polarization  antenna,  with 
vertical  and  horizontal  antennas  located  in  the  z  and  y  directions,  is  expressed  as 
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^(0=[^'l(0,>!*k0]r=E[4''£,(0-i4*1£,(0^r=2;[^;y-,(0,Wl(0cos«)]r, 

i=l  ;=1 


(12) 


where  E,(t)  as  the  electric  field  corresponding  to  the  ith  signal,  a]j]  and  a1/1 ,  respectively. 


are  the  /  th  elements  of  the  vertically  and  horizontally  polarized  array  vectors,  aM(^()  and 


aw  ($ ) ,  and  “  •  ”  represents  the  dot  product.  It  is  assumed  that  the  array  has  been  calibrated 


and  both  alvJ(^,)  and  a 1  (#,)  are  known  and  normalized  such  that 

|aM(4f  =|hW«)f  =  m  .  It  is  noted  that,  since  the  cos(^)  term  is  common  in  all  the 

elements  of  the  horizontally  polarized  array  manifold,  it  can  be  absorbed  in  the  array 
calibration  for  the  region  of  interest  and,  therefore,  removed  from  further  consideration. 
Then,  the  above  equation  is  simplified  to 

1,(0  =  [a^s{;\t),a^hf\t)}T  =  ,(0kurF  °  k-i,cf2 r )=  s,(t)Aa  o c, ,  (13) 


where  the  vector 

c,.=[c,,cJ=[cos(7/),sin(,,)^f  (14) 

represents  the  polarization  signature  of  the  i  th  signal. 

3.2  Polarlmetric  Time-Frequency  Distributions  (PTFDs) 

For  the  /  th  dual-polarized  sensor,  we  define  the  self-  and  cross-polarized  TFDs, 
respectively,  as 


Dxinx  5*.  (Of)  =  \ty(t  -  u,  r)x}'1  (u + (u  -  -))"  e~j2,!fr  dud  r  (15) 
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and 


-D  (.1^1  (t,f)  =  ~  u,T)x\i](u  +  ^)(x]k](u  -  r-^j)  e31* dudx ,  (16) 

where  the  superscripts  /  and  k  denote  either  v  or  h.  Define 

x,o) =[*,",,*rif ,  (H) 

to  include  the  two  polarization  components  received  in  a  dual-polarized  antenna.  Then,  the 
self-  and  cross-polarized  TFDs  constitute  the  2x2  polarimetric  TFD  (PTFD)  matrix 

D  (/,/)  =  f  \<p{t  -  u,  t)jc,  (u  +  — )xf  (u  -  —)e~j2!tfr dudr .  (18) 

2/J;  JJ  2  2 

The  diagonal  entries  of  D  (/,/)  are  the  self-polarized  TFDs  D  fj]  [M  (t,f) ,  whereas  the  off- 

JC  i  X  i  Xj  x{ 

diagonal  elements  are  the  cross-polarized  terms  D  ( t,f),i  &  k . 

X,  X, 

3.3  Spatial  Polarimetric  Time-Frequency  Distributions  (SPTFDs) 

Equations  (12)  -  (18)  correspond  to  the  case  of  a  single  dual-polarization  sensor.  With  an 
array  of  m  dual-polarized  antennas,  the  data  vector,  for  each  polarization  i ,  i  =  v  or  h ,  is 
expressed  as 

x['1(0  =  [x['1(0,x''1(/),...,4;!(0]r  =y[/](0  +  n[/](0  =  +  (19) 

Now,  we  generalize  the  PTFDs  to  a  multi-sensor  receiver.  Based  on  Eq.  (19),  the  following 
extended  data  vector  can  be  constructed  for  both  polarizations, 
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It  is  clear  that  the  dual-polarization  array,  compared  to  single-polarization  case, 
doubles  the  vector  space  dimensionality.  Specifically,  when  a[vj  (<fk )  =  a[/l1  {<j>k )  =  a(^, ) ,  the 
above  equation  simplifies  to 

»<&)  =  »W(A)®V  (27) 

It  is  now  possible  to  combine  the  polarimetric,  spatial,  and  t-f  properties  of  the  source  signals 
incident  on  the  receiver  array.  The  STFD  of  the  dual-polarization  data  vector,  x(?) ,  can  be 
used  in  formulating  the  following  SPTFD  matrix 

Dxx(t,/)  =  j]V(/  -  n,  t)\(u  +  ~)\H  (u  -  ^)e'jl¥r  dudt .  (28) 

This  SPTFD  matrix  embodies  the  information  required  to  address  typical  problems  in  array 
processing,  including  direction  finding,  as  will  be  shown  in  the  next  section. 

When  the  effect  of  noise  is  ignored,  the  SPTFD  matrix  is  related  to  the  source  TFD 
matrix  by  the  following  equation, 

Dxx(f,/)  =  B(0)QDss(f,/)Q"BH(0).  (29) 

4.  Polarimetric  Time-Frequency  MUSIC  (PTF-MUSIC) 

Time-frequency  MUSIC  (TF-MUSIC)  has  been  recently  introduced  for  improved  spatial 
resolution  for  signals  with  clear  t-f  signatures  [7].  It  is  based  on  the  eigenstructure  of  the 
STFD  matrix.  The  PTF-MUSIC  is  an  important  generalization  of  the  TF-MUSIC  to  deal 
with  diversely  polarized  signals  and  polarized  arrays.  It  is  based  on  the  SPTFD  matrix  of  Eq. 
(28). 

Consider  the  following  spatial  signature  matrix 
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(30) 


m= 


1  fa^) 


0 

•w(#). 


corresponding  to  DOA  <f> .  Since  al,|(^)j|  -m,  FH  (</>)¥  (0)  is  the  2x2  identity  matrix. 


To  search  in  the  joint  spatial  and  polarimetric  domains,  we  define  the  following 
spatio-polarimetric  search  vector 


w’c)=f§=FWc' 


(31) 


where  the  vector  c  =  [c,,c2]r  is  a  unit  norm  vector  with  unknown  polarization  coefficients 


c,  and  c2 .  In  Eq.  (31),  we  have  used  the  fact  that  ||F(^)c||  =  [c^F77  (^)F(^)cj2  =  (c/7c)2  =  1 . 


The  PTF-MUSIC  spectrum  is  given  by  the  following  function  [22,  23] 

/>(^)  =  [minf77(^,c)U„Uff(^,c)]-1  =[mincwFw(^)U„UfFWc]-1,  (32) 

c  c 

where  U„  is  the  noise  subspace  obtained  from  the  SPTFD  matrix  using  the  selected  t-f 
points.  Time-frequency  averaging  and  joint  block-diagonalization  are  two  known  techniques 
that  can  be  used  to  integrate  the  different  STFD  or  SPTFD  matrices  constructed  from 
multiple  t-f  points  [5,  7,  25].  The  selection  of  those  points  from  high  energy  concentration 
regions  pertaining  to  all  or  some  of  the  sources  enhances  the  SNR  and  allows  the  t-f  based 
MUSIC  algorithms  to  be  more  robust  to  noise  compared  to  its  conventional  MUSIC 
counterpart  [5], 

In  Eq.  (32),  the  term  in  brackets  is  minimized  by  finding  the  minimum  eigenvalue  of 
the  2x2  matrix  F//(^)UnU^F(^)  .  Thus,  a  computationally  expensive  search  in  the 
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polarization  domain  is  avoided  by  performing  a  simple  eigen-decomposition  on  the  2x2 
matrix.  As  a  result,  the  PTF-MUSIC  spectrum  can  be  expressed  as 

W  =  ^L[Fir(#)UIIU"F(^)],  (33) 

where  2min  [•]  denotes  the  minimum  eigenvalue  operator.  The  DO  As  of  the  sources  are 
estimated  as  the  locations  of  the  highest  peaks  in  the  PTF-MUSIC  spectrum.  For  each  angle 
(j)k  corresponding  to  one  of  the  n  signal  arrivals,  k  =  1,2,.  ..,w ,  the  polarization  parameters  of 
the  respective  source  signal  can  be  estimated  from 

KM  =  vmin[F"  (4)U„U„"F(^)],  (34) 

where  vmin[-]  is  the  eigenvector  corresponding  to  the  minimum  eigenvalue  2min  [•] . 


5.  Spatio-Polarimetric  Correlations 

In  a  conventional  single-polarized  array  problem,  the  spatial  resolution  capability  of  an 
array  highly  depends  on  the  correlation  between  the  propagation  signatures  of  the  source 
arrivals  [5,  26].  This  is  determined  by  the  normalized  inner  product  of  the  respective  array 
manifold  vectors.  In  the  underlying  problem,  in  which  both  the  spatial  and  polarimetric 
dimensions  are  involved,  the  joint  spatio-polarimetric  correlation  coefficient  between  sources 
/  and  k  is  defined  using  the  extended  array  manifold  a(^) ,  as 


P,t  )»(A  )  =  (aw  (A  if  »w(*  ) + ck2cn  (a1*1  ^  )f  a™  (*  )) 

m  m  ' 


(35) 


—  Ck\Cl\P),k  +  Ck2CnP\jk 
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where 


m 


(</>/)  is  the  spatial  correlation  coefficient  between  sources  /  and 


k  for  polarization  i ,  with  i  denoting  v  or  h. 


An  interesting  case  arises  when  the  vertically  and  horizontally  polarized  array 
manifolds  are  identical,  i.e.,  aw(^)  =  a In  this  case,  p)'^  =  p)hk ,  and  the  joint  spatio- 

polarimetric  correlation  coefficient  becomes  the  product  of  the  individual  spatial  and 
polarimetric  correlations,  that  is, 


Pit  =  Pit  Pi, 


with 


Pijc  =  cf c/  =  cos(y/)cos(yA)e7('?'“%)  +sin(y/)sin(yi)  (37) 

representing  the  polarimetric  correlation  coefficient.  In  particular,  for  linear  polarizations, 
tj /  =fjk= 0 ,  and  Eq.  (37)  reduces  to 


p,k=cos(y,-yk).  (38) 

Since  |p;  k  <  1 ,  with  the  equality  holds  only  when  the  two  sources  have  identical  polarization 

states,  the  spatio-polarization  correlation  coefficient  is  always  equal  to  or  smaller  than  that  of 
the  individual  spatial  correlation  coefficient.  The  reduction  in  the  correlation  value  due  to 
polarization  diversity,  through  the  introduction  of  pl  k ,  translates  to  improved  distinctions  of 
sources.  As  such,  two  sources  that  could  be  difficult  to  resolve  using  the  single-polarized 
spatial  array  manifold  a[v](^)  or  a[/,](^)  can  be  easily  separated  using  the  extended  spatio- 
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polarized  array  manifold,  defined  by  a(^) .  This  improvement  is  more  evident  in  the  case 
when  the  source  spatial  correlation  is  high,  but  the  respective  polarimetric  correlation  is  low. 


6.  Moving  Sources  with  Time-Varying  Polarizations 
In  this  section,  we  consider  the  performance  of  tracking  moving  targets  with  time- varying 
polarization  signatures.  Time- varying  polarizations  are  often  observed  when  active  or  passive 
sources  move  or  change  orientations  [27].  The  performance  of  polarimetric  MUSIC  (P- 
MUSIC)  and  PTF-MUSIC  techniques  are  discussed  and  compared.  We  show  that,  the 
SPTFD  maintains  the  instantaneous  polarization  information  at  the  auto-term  t-f  points, 
whereas  the  energy  of  desired  source  signals  is  localized.  Careful  treatment  of  the  data 
samples  is,  however,  required  to  effectively  utilize  the  polarization  diversity,  since  the 
construction  of  the  covariance  and  SPTFD  matrices  using  all  available  data  samples  may 
compromise  the  polarization  distinctions  between  sources. 

Because  the  SNR  enhancement  using  multi-sensor  t-f  processing  is  well  documented 
[5,  6,  9],  we  consider  the  noise-free  case  in  this  section  for  notational  simplicity.  In  this  case, 
the  received  signal  vector  is 


*(0  = 


A",(®(0XqI"(0°aM)" 

A,,|(<&(0)(q,‘l(0°s(0). 


(39) 


with  A">(<I>(0)  =  [a|,|W1(0),  ....  *“'(«>„(')>] =  v,h  is  the  time- varying  array  response 
matrix.  Note  that  we  used  <X>(/)  =  . . . ,  to  emphasize  the  fact  that  the  DO  As  are 

now  time-varying. 
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6.1  Polarization  Diversity  and  Selection  of  Data  Samples 

In  order  to  properly  utilize  polarization  diversity,  careful  attention  should  be  paid  in  the 
selection  of  data  samples.  We  use  the  P-MUSIC  as  the  example. 

Given  the  time- varying  nature  of  the  source  signal  polarizations,  the  covariance  matrix  of 
the  received  signal  vector  is 


Rxx(t)  =  E[x(t)xH(t)]  =  E 


A1*1  («(f  )Xq,*'(0  °  *('))!  a1*1  (®(OXqw  (')  °  *>(<)). 


(40) 


When  the  covariance  matrix  is  replaced  by  using  a  time  average  over  the  entire  data  samples 
and  considering  the  special  case  where  the  spatial,  temporal,  and  polarimetric  signatures  can 
be  decoupled,  it  becomes 


qw(')(q'"(t))"  °R, 
<T«)(qw(0)*  °R» 


ql"(tXg'“(0)1'  R 
ql'"(0(ql‘,«))"  °R 


ql"(0(q"'(0)"  ql’l(0(q'*l(0)" 

Rss 

Rss  ] 

.q'*l(')(q,’l(0)"  q'*'('Xqw (<))". 

Rss 

Rss  j 

(41) 


where  Rss  =  s(t)sH  (t) ,  and  B  = 


A[vl($(0)  0 

0  A1*1  (0(0) 


The  time- varying  source  signal 


polarization  vectors  are  defined,  similarly  to  Eqs.  (23)  and  (24),  as 


qlv]  (0  =  [cos(y,  (0),  •  -  - ,  cos(/„  (t))]7 


q[A1(t)  =  [sin(y,  {t))eJnM , . .  .,sin(y„  (t))ejnM  f . 
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Consider  an  alternative  scenario  where  the  polarimetric  signature  is  time-invariant.  Then  the 
covariance  matrix  can  be  expressed  as 


qw(qM)" 

qw(qw)" 


q[v](qtA])" 

qw(qw)tf 


o 

Rss 

Rss 

Rss 

Rss 

(42) 


The  effect  of  the  signal  time- varying  polarization  on  the  covariance  matrix  is  clear  from  Eqs. 
(41)  and  (42).  When  q[,I(0(qtAI(0)ff  -ql,!(qW)H  for  i,k  =  v,h  ,  the  corresponding 
covariance  matrices  are  identical,  and  the  two  cases  of  sources  with  time-varying  and  time- 
invariant  polarizations  will  lead  to  the  same  performance.  Consider,  for  example,  two 
linearly  polarized  sources  whose  polarization  angle  change  linearly  over  [0°,90°]  for  yx  and 
[90°,0°]  for  y2 .  The  DOAs  are  assumed  time-invariant.  This  case  is  equivalent  to  both 
sources  assuming  fixed,  time-invariant  polarization  of  y  =  45  °,  and  thereby,  the  source 
polarization  diversity  cannot  be  utilized  in  DOA  estimation  using  P-MUSIC. 

To  utilize  the  source  signal  polarization  diversity,  the  received  data  covariance  matrix 
in  (41)  should  be  constructed  using  a  selected  set  of  data  samples  (e.g.,  using  moving  average 
of  the  received  data  samples)  such  that  the  polarization  distinction  between  sources  is 
maintained,  even  if  the  source  DOAs  show  insignificant  changes  over  the  observation  period. 
In  general,  for  moving  target  tracking  problems,  the  selection  of  the  data  window  for  the 
moving  average  covariance  matrix  should  be  determined  such  that  the  DOA  variation  over 
the  window  is  small. 

We  maintain,  however,  that  the  reduction  of  data  samples  may  compromise  the 

subspace  estimation  and,  subsequently,  DOA  estimation  performance,  particularly  when  the 

SNR  is  relatively  low,  and/or  the  sources  are  closely  spaced.  The  use  of  t-f  representations 
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overcomes  this  shortcoming  and  provides  robust  subspace  estimation  in  rapidly  time-varying 
polarization  environments. 

6.2  Polarization  Diversity  Consideration  in  PTF-MUSIC 

We  now  show  that  the  instantaneous  polarization  information  is  incorporated  in  the 
SPTFD  and,  therefore,  PTF-MUSIC  maintains  polarization  diversity,  while  providing  SNR 
enhancement  and  source  discriminations. 

In  the  presence  of  time-varying  polarized  sources,  the  auto-  and  cross-polarized  SPTFD 
matrix  can  be  expressed  as 


(q[f,(«+-)(qW(w--)f)°(s(«+-)(s(» 


e  duck 


=  Jj$/ - u, rJA1'1  (<$(t + -^[g1*1  (m, r) ° K(«, r)^A[*] (<5(t - e™Tdudc 


(43) 


where  G['k]  ( t ,  r)  =  q[,]  (t  +  -^Xq1^1  (t  -  ^j)H ,  and  K(t,  r)  =  s(t  +  ~)sn  (t  -  .  When  the  window 

length  of  the  applied  t-f  is  small, 


(44) 

where 

~  u,  r)G[,*]  (u,  t)  °  K(w,  r)e^j27fT  dndr .  (45) 


We  assume  that  the  frequency  and  the  polarization  signatures  of  the  sources  change  almost 

linearly  within  the  temporal  span  of  the  t-f  kernel.  Then,  using  the  first-order  Taylor  series 
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expansion,  the  polarization-dependent  terms  can  be  approximated  as 
XT  d 

Vi(t +—)  =  yi(t)+—yi( t),  where  y,(t)  = — The  auto-terms  of  the  source  polarization 


information,  which  reside  on  the  diagonals  of  G1”1'1  (/,r)  ,  G[W,](t,r)  ,  G[Av](t,r)  and 
G[hh](t, t),  are  given  by 


[G[vv)(/,r)],  =  ^[cos(2/,.(0)  +  cos(z7,(/))] 

(46) 

[g[v/,1(/,t)1,  =  ^[sin(2//(0)-sin(T7i(0)] 

(47) 

[gc*v,(/,t)1(  =  ^[sin(2y,(0)  +  sin(T7,(0)J 

(48) 

?rlhh](t,T)\  =  ^(-cos(2y,(0)  +  cos(T7,(/))J , 

(49) 

respectively.  The  second  sinusoidal  terms  in  Eqs.  (47)  and  (48)  results  in  zero  values  in 
Dpv for  symmetric  t-f  kernels.  Therefore,  Ds(fy*](*,/)can  be  expressed  at  the  auto¬ 
term  points  as 


with 


Djp#\  (t’f) = |cos(2  rM°SiSi  (!,f) + c, 7(/,/) 


(50) 


=  -^cos(2y,(0)^j(  (t,f)  +  ca(t,f) 


(51) 


( t,f )  =  Dpp  0 t,f )  =  —  sin(2/<(t))£)  ^  (/,/) 


(52) 
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cu{t,f)  =  |  - u, t ) cos(z7, (/ ))[K(7 , r )](7 e  j24Tdudr . 


(53) 


When  different  sources  are  uncorrelated,  their  t-f  signatures  have  no  significant  overlap.  If 
the  t-f  points  located  in  the  auto-term  region  of  the  i  th  source  are  used  in  constructing  the 
SPTFD  matrix,  then 


D  „(/,/)  = 


»M(f*,(0)  0 

0  a™(fi,(/)) 


M. 


a'"«(0)  0 

0  a'*'(«S,(0) 


H 


(54) 


where 


M/  =  — D  (t,f) 

cos(2y,(t))  sin(2y,(0)  " 

+ 

1 

JTJ 

o 

_ 1 

sin(2y,(0)  -cos(2y,.(0). 

[ 

o 

In  the  new  structure  of  the  SPTFD  matrix  of  Eq.  (54),  the  source  time-varying  polarization 
has  the  effect  of  loading  the  diagonal  elements  with  c„(t,f)  and,  as  such,  alters  the 

eigenvalues  of  the  above  2x2  matrix.  However,  the  eigenvectors  of  My  remain  unchanged. 

The  new  eigenvalues  are  Xi  l  2  —  cn  {t,  f)  ±  ■—  Dss  (/,  /) .  The  signal  polarization  signature,  i.e., 

the  eigenvector  corresponding  to  the  maximum  eigenvalue,  is 
v/>ax  =  [cos(y,(/)),  sin(y/(t))f .  Therefore,  in  the  context  of  PTF-MUSIC,  the  instantaneous 
polarization  characteristics  can  be  utilized  for  source  discriminations. 


7.  Simulation  Results 

Two  chirp  signals  are  considered  and  their  parameters  are  listed  in  Table  I.  The  two 
sources  impinge  on  a  uniform  linear  array  (ULA)  of  five  cross-dipoles  with  an  interelement 
spacing  of  half  a  wave-length.  The  array  responses  in  both  horizontal  and  vertical 
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polarizations  are  identical.  The  input  SNR  is  5dB.  The  source  signals'  polarization  angles 
(Vi(t),i  =  1,2)  change  linearly  with  time  over  the  observation  period  of  T  =  256  samples 

(Figure  2).  Note  that  sources  with  polarization  angles  of  0°  and  180°  are  both  horizontally 
polarized,  whereas  a  polarization  angle  difference  of  90°  implies  orthogonal  polarizations. 
Figure  3  shows  the  polarization-averaged  pseudo  Wigner-Ville  distribution  (PWVD)  of  the 
received  data  at  the  reference  sensor  with  a  window  length  of  65.  When  the  sources  have 
low  polarization  correlation,  the  cross-terms  of  the  signal  components  corresponding  to  the 
vertical  and  horizontal  polarizations  cancel  out,  resulting  in  low  polarization-averaged  cross¬ 
term  values. 

In  tracking  the  moving  sources'  DOAs,  different  sets  of  SPTFD  matrices  are 
constructed.  Each  set  uses  P  consecutive  (neighboring)  t-f  points  of  the  two  sources'  auto¬ 
terms,  with  the  middle  t-f  point  being  at  tt  .  The  objective  is  to  examine  the  proposed 

algorithm  performance  at  different  source  polarization  states  and  to  demonstrate  the  tracking 
accuracy  of  the  algorithm.  This  is  achieved  by  applying  the  PTF-MUSIC  for  each  set  of 
SPTFD  matrices  separately. 

Figure  4(a)  shows  the  tracking  accuracy  of  the  PTF-MUSIC  algorithm  for  the  DOA 
estimations  of  the  two  sources  with  P-3 1  and  using  a  PWVD  with  window  length  1=65.  The 
performance  of  P-MUSIC  is  also  shown  in  Figure  4(b)  for  P= 31  and  in  Figure  4(c)  for  P-95. 
The  tracking  accuracy  of  the  algorithms  was  evaluated  at  49  time  points  that  are  spaced  5 
samples  apart  (i.e.,  6=8  and  t49  =  248 ).  We  note  that  both  techniques  use  the  same  data 

samples.  The  PTF-MUSIC  generally  outperforms  the  P-MUSIC  as  can  be  seen  in  Figure  4. 
The  performance  of  the  PTF-MUSIC  is  particularly  impressive  in  regions  where  the  signal 
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energy  is  highly  localized  in  the  t-f  domain  and  the  polarization  distinction  between  the  two 
sources  is  high.  The  importance  of  polarization  diversity  is  emphasized  by  the  fact  that,  in  the 
central  and  edge  regions  where  the  two  sources  have  low  polarization  distinction,  the 
performance  of  both  PTF-MUSIC  and  P-MUSIC  suffers. 

8.  Conclusion 

This  paper  considered  the  application  of  polarimetric  time-frequency  (PTF-MUSIC) 
for  tracking  moving  sources  with  time-varying  polarizations.  It  has  been  shown  that  the 
difference  in  the  instantaneous  polarizations  of  the  sources  can  be  uniquely  utilized  by  the 
proposed  approach  to  maintain  polarization  diversity.  It  is  shown  that  the  PTF-MUSIC 
approach  is  superior  to  a  covariance  matrix  approach  that  is  based  on  a  moving  window  of 
the  received  data.  The  proposed  approach  also  allows  for  the  discrimination  of  sources  based 
on  their  time-frequency  signal  characteristics. 
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Table  I.  Signal  parameters  of  moving  sources  with  time- varying  polarizations 


Start 

End 

DOA 

7 

n 

freq. 

Freq. 

(deg-) 

(deg.) 

(deg.) 

Source  1 

0 

0.20 

-15  to  5 

Oto  180 

0 

Source  2 

0.20 

0.40 

-10  to  10 

180  to  0 

0 
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Time  (samples) 


(a)  PTF-MUSIC  (using  31  time-frequency  points) 


Time  (samples) 


(b)  P-MUSIC  (using  31  data  samples) 


Time  (samples) 

(c)  P-MUSIC  (using  95  data  samples) 

Figure  4.  PTF-MUSIC  and  P-MUSIC  tracking  performance. 
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DOA  and  Polarization  Estimation  for 
Wideband  Sources 

Baha  A.  Obeidat,  Yimin  Zhang,  and  Moeness  G.  Amin 

Abstract 

This  paper  presents  an  approach  to  direction-of-arrival  (DOA)  and  polarization  estima¬ 
tions  of  far-field  wideband  signals  using  an  array  of  cross-polarized  dual-feed  sensors.  The 
wideband  signals  are  decomposed  into  a  set  of  narrowband  signals.  Coherent  signal  sub¬ 
space  processing  is  then  applied  to  dimension-extended  source  signal  correlation  matrices 
at  different  narrowband  frequencies,  leading  to  the  construction  of  single  signal  subspace. 
The  utilization  of  polarization  diversity  through  the  use  of  a  dual-feed  cross-polarized  ar¬ 
ray  to  process  the  polarized  signal  arrivals  lends  an  additional  degree  of  freedom,  leading 
to  performance  improvement  over  the  case  where  a  single-polarization  array  is  used. 
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I.  Introduction 


Direction-of-arrival  (DOA)  estimation  for  far-field  sources  has  been  the  subject  of  nu¬ 
merous  array  processing  studies  [1]  [2],  When  the  signals  that  are  incident  on  an  array  are 
wideband,  coherent  signal  subspace  processing  [3]  can  be  used  to  apply  typically  narrow- 
band  subspace-based  estimation  methods  to  the  wideband  case.  Coherent  signal  subspace 
processing  relies  on  the  frequency  domain  decomposition  of  the  array  output  into  a  number 
of  narrowband  outputs  to  which  the  conventional  narrowband  array  processing  model  ap¬ 
plies.  For  each  frequency  subband,  transformation  matrices  are  designed  that  focus  their 
respective  covariance  matrices  to  a  particular  frequency  of  interest.  The  covariance  ma¬ 
trices  are  later  coherently  combined  to  produce  a  coherent  signal  space  [4].  This  approach 
outperforms  its  non-coherent  counterpart. 

Polarization  and  polarization  diversity,  on  the  other  hand,  have  been  used  in  wireless 
communication  systems  [5]  [6]  and  in  various  types  of  radar  systems  [7].  Additionally,  they 
have  been  incorporated  into  array  processing  for  DOA  estimation  for  narrowband  sources 
[8].  The  incorporation  of  signal  polarization  information  lends  an  additional  degree  of 
freedom  and  consequently  results  in  improved  DOA  estimation  and  signal  synthesis  [9] 
[10]. 

This  paper  considers  the  estimation  of  DOA  and  polarization  signatures  of  wideband 
sources.  We  show  that  the  coherent  subspace  processing  remains  effective  in  dual-polariza¬ 
tion  applications.  Because  the  array  manifold  in  one  polarization  may  differ  to  that  of  the 
other  polarization,  the  transformation  matrices  in  general  should  be  design  separately  for 
different  polarizations.  The  use  of  polarization  diversity  is  particularly  effective  when  the 
source  signals  are  closely  spaced  but  have  distinct  polarizations.  Simulation  results  are 
provided  which  show  the  superior  performance  of  the  proposed  method  in  comparison  to 
single-polarized  arrays  with  identical  as  well  as  extended  apertures. 

II.  Signal  Model 

A.  Signal  Model 

Consider  an  M-element  array  of  identical  cross-polarized  (such  as  horizontal  and  verti¬ 
cal),  dual-feed  sensors  as  shown  in  Fig.  1.  The  vertically  polarized  sensors  are  placed  in 
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the  z-axis  and  the  horizontally  polarized  sensors  are  placed  in  the  y-axis.  Upon  the  array 
N  polarized  farfield  wideband  sources  are  incident.  It  is  assumed  that  the  source  signals 
lie  in  the  x-y  plane  such  that  their  elevation  angles,  6n  are  90  degrees  and  their  azimuth 
angles  are  denoted  as  4>n,  n  =  1, 2, . . . ,  N.  The  source  signals  are  assumed  to  have  fixed 
polarizations  over  the  observation  period.  Additionally,  since  the  sources  are  farfield,  it  is 
practical  to  assume  that  the  source  signal  polarizations  do  not  change  as  received  by  one 
sensor  to  the  next. 

The  received  data  across  the  M  vertically  polarized  sensors  is  converted  to  baseband, 
time-sampled  as 

xM(t)  =  AM(i)sM(i)+nM(t)  (1) 

with  the  array  response  matrix  as 

AM(*)  =  aW (</>i, t)  ...  a [v](<f>N,t)  (2) 

where  a^l (<pn,t)  denotes  the  array  response  vector  for  the  nth  source  signal,  n^l (t)  is  the 
temporal  domain  additive  sensor  noise  vector.  The  received  signal  vector  in  the  vertical 
polarization  is  s^l  (t)  which  can  be  decomposed  as  sM(t)  =  Q^s(t)  since  the  source  po¬ 
larization  is  constant  over  the  observation  period  and  independent  of  the  frequency,  s  (t) 
is  now  the  polarization-free  received  signal  vector  and  the  matrix  containing  the  source 
signal  vertical  polarization  components  as 

QH  =  diag[cos(7i),  . . .  ,  cosfrw)].  (3) 

Then,  (1)  can  be  rewritten  as 

XH  (t)  =  AM  (t)QMs(t)  +  nW  (t)  (4) 

The  output  of  the  array  is  then  frequency  decomposed  into  J  nonoverlapping  complex 
subband  components.  The  bandwidth  of  each  subband  is  much  smaller  than  its  central 
frequency  fj,j  =  1,2 , ...,  J,  so  that  each  subband  snapshot  approximately  satisfies  the 
classical  narrowband  assumption  [4].  The  kth,  k  =  1, 2, . . . ,  K,  vertically  polarized  snap¬ 
shot  received  signal  vector  follows  the  narrowband  model  and  is  written  as 

4 fUi)  =  AWt/jJQHstt/j) +  »«(/,)  (5) 
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where 


(6) 


AM(/j)  =  [a ■■■, *n  (<Pn,  fj)\ 

is  the  vertical  array  response  matrix  for  the  frequency  fj  and  polarization  with  a n\<t>n,  fj) 
as  the  vertical  array  response  of  the  nth  source  at  the  jth  frequency,  nj %\fj)  is  the  Mxl 
noise  vector  in  the  vertical  polarization  at  kth  sample  of  the  frequency  fj.  s k(fj)  is  the 
polarization-free  source  signal  vector  at  frequency  fj. 

The  same  type  of  manipulation  can  be  applied  to  the  received  signal  across  the  hori¬ 
zontally  polarized  array  to  yield 

=  Alh](fj)Q.[h]sk(fj)  +  nf](fj),  (7) 

with 

=  diag[sin(71)eJ,n,  ...  , sin(7Ar)ejr,N]  (8) 

Concatenating  the  received  signal  in  the  horizontal  and  vertical  polarizations  provides  the 
extended-dimension  received  signal  vector, 

xtHfi)]  =  k*(/,)l 

xfOi)]  [AWyjJQWstl/,)]  [4A'(/i)J  ' 

III.  Proposed  Method 

The  idea  of  coherent  signal  subspace  was  first  proposed  in  [3].  When  A W  (/,-),  i  —  v  or 
h,  j  =  1, 2, . . . ,  J,  are  full  rank,  one  can  design  the  transformation  matrices,  T such 
that 

T'i|(/i)Al,(/j)  =  AB(/0).  (10) 

In  the  underlying  dual-polarization  scenario,  it  can  be  cast  in  a  block  diagonal  form  to 
give 

T(/ )  fAM(/j)l  A  fTM(A)  0  1  [am(/j)1  =  [aW(/o)1  (11) 

LAWte>J  l  0  TW(t)J  [aW(/,)J  [aW(/o)J‘ 

Generally,  since  the  array  responses  in  the  vertical  and  horizontal  polarization  may  be 
different,  the  two  focusing  matrices  are  designed  to  focus  the  array  response  matrix  in  the 
vertical  or  horizontal  polarizations  to  the  frequency  of  interest. 
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Applying  the  block  diagonal  transformation  matrix  to  the  extended-dimension  received 
signal  vector  will  yield 

y(fj)  =  n/Mfi)  (is) 

from  which  a  new  covariance  matrix  can  be  constructed 

j 


Ryy  =  (T(fMfi)) 

3= 1 

AM(/0)qw1  .j 

AW(/o)QW 


H 


3=1 


AM(/o)QM 

AW(/o)QW 


H 


+  a2J^T(fi)P(fj)THUi) 

j= 1 


AM(/q)QM 

AW(/o)QW 


RL 


AM(/o)QM 

AW(/0)qW 


j? 


+  <7^P 


2-rV 


where 


and 


RL  =  £r»(/3) 

3= 1 


p'=Ew>)p(/i)T,,(/j) 

i=i 


(13) 


(14) 


(15) 


A.  Spatio-polarimetric  correlations 


It  is  evident  from  (13)  that  the  array  has  a  extended-dimension  aperture.  In  this  struc¬ 
ture,  the  closeness  of  two  sources,  l  and  m,  is  determined  by  their  spatio-polarimetric 
correlation  coefficient.  Consider  a  simple  scenario  where  the  array  responses  at  the  focus¬ 
ing  frequency  are  identical  in  both  polarizations,  i.e.  a^(0, /0)  =  a^(0, /0)  =  a(</>,  /0), 
the  spatio-polarimetric  correlation  coefficient  is  defined  as  the  product  of  the  spatial  cor¬ 
relation  coefficient  and  the  polarimetric  correlation  coefficient  [11],  that  is 

An,i  —  Pm,lPl,m •  (16) 

In  (21),  the  spatial  correlation  coefficient  is  the  normalized  inner  product  of  the  sources’ 
array  responses  at  the  focusing  frequency, 


M 


a£(&/o)aj(&/o), 
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and  their  respective  polarimetric  correlation  coefficient  is 

PmjL  =  cos(7;)cos(7m)ej(,!l_7!m)  +  sin(7()sin(7m).  (18) 

In  particular  when  the  soruces  are  linearly  polarized,  rjn  =  0,  n  —  1, 2, . . . ,  iV,  (19)  is 
simplified  as 

PmjL  =  cos(7i  ±  jm).  (19) 

When  an  M-e lement  array  of  non-polarization  sensitive  sensors  is  used,  the  spatial  cor¬ 
relation  coefficient  between  two  source  incident  upon  it  is  merely  j3mti.  Since  \pmji\  <  1, 
the  use  of  a  dual-polarized  array  can  only  reduce  the  spatio-polarimetric  correlation  coef¬ 
ficient.  The  reduction  is  particularly  significant  when  the  sources  are  closely  spaces  but 
have  distinct  polarizations. 

A  2M-element  array  of  equally  spaced  non-polarization  sensitive  sensors  will  have  a 
spatial  correlation  of 

Pm, l  "  Pm, l Pi, mi  (20) 

where  pitTn  is  the  reduction  is  in  the  spatial  correlation  caused  by  the  doubling  the  number 
of  sensors  and  extending  the  array  aperture.  This  coefficient  is 

d 

|/S;,m|  =  cos  7T M  —  (sin (j>i  sin  (f>m )  ,  (21) 

where  Aq  is  the  wavelength  of  the  focusing  frequency. 


B.  Polarimetric  MUSIC 

The  coherent  signal  subspace,  Es  =  je1;  ...  ,  e^j ,  of  dimension  2 M  x  N  is  con¬ 
structed  from  the  eigenvectors  corresponding  to  the  N  largest  eigenvalues  of  Ryy.  The 
noise  subspace,  on  the  other  hand,  is  constructed  using  the  smallest  2M—N  and  is  denoted 
En.  The  polarimetric  MUSIC  spectrum  is  given  by 


p«  = 


where  Amjn  denotes  the  minimum  eigenvalue  operator  and 


B(0,/o) 


aw(^/o)  0 

0  a^(<^,/o) 


(23) 
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is  the  dual-polarization  search  matrix  at  the  focusing  frequency  /0. 

The  DOAs  of  the  sources  are  estimated  as  the  locations  of  the  highest  peaks  in  the 
MUSIC  spectrum.  For  each  angle  (j)n ,  n  —  1, . . . ,  N,  the  polarization  parameters  of  the 
respective  source  signal  can  be  estimated  from 

=  vmin[BH(^,/„)B„E"B ■(&/„)],  (24) 

where  vmin[-]  is  the  eigenvector  corresponding  to  the  minimum  eigenvalue  Amjn[-]. 

IV.  Simulations 

A  four-sensor  uniform  linear  cross-polarized  array  is  considered  with  each  sensor  con¬ 
sisting  of  a  pair  of  horizontally  and  vertically  polarized  elements.  The  array  has  identical 
responses  in  the  horizontal  and  vertical  directions.  The  interelement  spacing  is  half  of 
the  wavelength  at  the  focusing  frequency,  /0  =  100  Hz.  Two  farfield  sources  with  DOAs 
of  (f>i  =  10°  and  </>2  =  11°  impinge  on  the  array  with  polarization  angles  71  =  45°  and 
72  =  —45°,  respectively.  The  two  sources  are  linearly  polarized,  i.e.,  their  polarizations 
phase  angles  ?h  =  r/2  =  0°.  The  two  sources  have  the  same  central  frequency,  /o  =  100  Hz, 
and  a  bandwidth  of  40  Hz.  The  SNR  is  set  to  OdB,  and  the  number  of  frequency  subbands 
used  was  J  =  33  with  each  of  K  =  100  snapshots.  The  search  grid  used  for  the  MUSIC 
and  polarimetric  MUSIC  methods  was  0.01°.  The  performance  of  the  polarimetric  MU¬ 
SIC  is  compared  to  that  of  the  conventional  MUSIC  using  a  similar  array  with  horizonal 
polarized  sensors  in  place  of  the  cross-polarized  sensors. 

Since  the  dual-polarized  requires  twice  the  processing  as  that  of  a  non-polarized  array 
with  4  sensors,  a  fairer  comparison  would  be  against  that  of  an  array  with  twice  the  number 
of  sensors,  i.e.,  8  sensors.  For  this  purpose  we  also  performed  simulations  to  examine  the 
performance  of  the  proposed  method  when  compared  to  that  of  a  non-polarized  array  and 
that  of  a  non-polarized  uniform  linear  array  with  twice  the  number  of  sensors.  It  is  worth 
noting,  however,  that  the  8-element  array  uses  more  than  twice  the  aperture  as  that  of 
the  dual-polarized  array  used  since  both  has  the  same  interelement  spacing.  In  fact,  the 
dual-polarized  array’s  aperture  is  while  the  8-element  array  has  an  aperture  of 

Figure  2  shows  a  single  trial  run  of  the  spatial  spectra.  Note  that  in  this  case  the 
two  sources  have  orthogonal  polarizations  and,  irrespective  of  the  spatial  separation  of 
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the  sources,  the  joint  spatio-polarimetric  correlation  coefficient  between  the  two  sources 
is  zero.  As  a  result,  while  the  conventional  MUSIC  cannot  resolve  the  two  closely  spaced 
sources,  the  polarimetric  MUSIC  provides  accurate  DOA  estimations. 

Figure  3  shows  the  root  mean  squared  error  (RMSE)  performance  of  the  DOA  estimation 
of  the  three  types  of  arrays.  The  results  are  obtained  from  100  independent  trials.  It  is 
evident  that  the  polarimetric  methods  outperforms  the  two  other  methods.  In  this  example 
the  spatial  correlation  coefficient  between  the  two  source  is  /?ii2  =  0.9982,  a  high  value 
due  to  their  close  spatial  separation.  However,  the  polarimetric  correlation  coefficient 
was  pit2  =  0  due  to  the  fact  that  the  two  sources  are  orthogonal  in  polarization.  That 
reduces  the  spatial  polarimetric  correlation  to  the  product  of  the  two,  i.e.,  zero.  As  for  the 
8-element  array,  the  spatial  correlation  it  its  case  is  /31)2  =  0.9924  which  remains  a  high 
value.  The  difference  in  the  spatial  correlation  or  the  spatial  polarimetric  correlation  can 
explain  the  difference  in  the  performance  of  the  three  methods. 

When  the  polarization  angle  of  source  1  is  —  30  and  that  of  source  2  is  71  =  60, 
whereas  all  other  variables  remain  unchanged,  the  performance  of  the  three  types  of  arrays 
is  compared  in  Fig.  4.  Here  the  polarimetric  correlation  coefficient  between  the  two  sources 
is  Pi,2  =  0.866  which  results  in  a  spatio-polarimetric  correlation  coefficient  of  0.8644  as 
opposed  to  0  in  the  previous  set.  The  performance  of  the  polarimetric  MUSIC  worsens  but 
outperforms  that  os  the  other  types  of  arrays.  As  in  other  polarimetric-based  methods,  the 
performance  of  this  method  in  DOA  estimation  is  a  function  of  the  polarization  diversity 
among  the  source  signals.  The  higher  the  polarization  diversity  the  better  the  performance. 

V.  Conclusion 

This  paper  proposes  a  direction-of-arrival  (DOA)  estimation  method  for  wideband  po¬ 
larized  farfield  sources  that  are  incident  on  an  array  of  cross-polarized  double-feed  sensors. 
This  method  utilizes  the  received  signal  across  the  array’s  two  polarization  arms  to  con¬ 
struct  a  dimension  extended  received  signal  vector.  Matrix  focusing  is  then  used  for  the 
purpose  of  constructing  a  coherent  signal  subspace.  When  the  source  signals  have  di¬ 
verse  polarizations  this  method  is  shown  to  outperform  the  conventional  non-polarimetric 
method. 
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Fig.  3.  RMSE  of  the  DOA  estimation. 


Fig.  4.  RMSE  of  the  DOA  estimation. 
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Sensor  Configuration  in  Polarized  Antenna 

Arrays 


Baha  A.  Obeidat,  Moeness  G.  Amin,  Yimin  Zhang,  and  Ahmad  Hoorfar 


Abstract 

Single-  and  dual-polarized  array  configurations  are  often  used  in  imaging  and  radar 
applications.  In  this  paper,  spatio-polarimetric  correlation  are  used  as  a  framework  to 
evaluate  the  direction  finding  performance  of  different  receiver  configurations  employing 
double-feed  dual-polarized  array  as  well  as  arrays  of  single-feed  single-polarized  elements. 
The  spatio-polarimetric  correlation  coefficients  are  derived  for  several  array  reconfigura¬ 
tions  and  used  to  provide  a  configuration-based  performance  comparison  based  on  spatial 
resolution  and  grating  lobes. 
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I.  Introduction 


It  is  well  known  that  the  use  of  dual  polarizations  can  enhance  the  performance  of 
direction-of-arrival  (DOA)  estimation  in  array  processing  [1],  [2],  [3]  and  improve  robust¬ 
ness  to  channel  fading  in  wireless  communications  [4]. 

Many  of  modern  communication  and  radar  systems  employ  antenna  arrays  made  of 
microstrip  or  other  printed  radiating  elements.  These  antenna  elements  have  many  ad¬ 
vantages,  including  low  profile,  lightweight,  easy  integration  to  various  shapes  and  surfaces 
and  low-production  cost.  One  main  disadvantage  of  microstrip  patch  antenna,  however, 
is  its  narrow  bandwidth.  Over  the  years,  many  techniques  have  been  reported  in  the  lit¬ 
erature  for  bandwidth  enhancement  of  these  low-profile  antennas  [5].  Extension  of  such 
wideband  techniques  to  design  of  dual-polarized  dual-feed  low-profile  antennas,  however,  is 
not  straightforward  if  one  also  requires  low  polarization  isolation  between  the  two  output 
(input)  ports,  which  is  often  the  case  in  many  sensitive  applications.  In  fact,  a  typi¬ 
cal  polarization  isolation  constraint  of  less  than  25dB  drastically  reduces  the  operating 
bandwidth.  One  possible  solution  to  avoid  the  polarization  isolation  issue,  and  the  as¬ 
sociated  complexity  of  two-port  elements  in  a  large  array  scenario,  is  to  use  a  wideband 
linearly  polarized  single-port  antenna  that  also  senses  the  orthogonal  polarization  by  a 
90°  mechanical  rotation  of  the  antenna  element  using  a  motor-driven  turntable.  Another 
alternative  solution  may  be  the  use  of  active  elements,  e.g.,  varactor  or  pin  diodes,  to 
control  and  switch  between  the  two  states  of  polarization  of  a  wideband  dual-polarized 
antenna  element  [6]. 

Therefore,  the  capability  to  reconfigure  a  double-feed  dual-polarized  array  into  a  sim¬ 
plified  array  of  single-feed  single-polarized  sensors  becomes  of  particular  interest  to  array 
processing  for  polarized  signals.  The  motivation  to  simplify  the  array  structure  in  hard¬ 
ware  stems  from  the  fact  that  single  polarized  antennas  and  antenna  arrays  are  more 
commercially  available  than  the  dual-polarized  double-feed  arrays.  Also,  with  the  capabil¬ 
ity  to  change  the  sensor  polarizations  within  the  array  it  becomes  important  to  compare 
the  performance  of  a  dual-polarized  array  with  its  single-polarized  sensor  reconfigurations. 

In  this  paper,  we  compare  the  performance  of  a  double-feed  dual-polarized  array  to 
other  affiliate  reconfigurations  consisting  of  single-feed  single-polarized  antenna  elements. 
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Of  interest  are  two  array  reconfigurations;  a  reconfigured  array  of  all  vertical  polarized 
antennas  which  is  concatenated  with  an  array  of  all  horizontal  polarizations,  and  a  recon¬ 
figured  array  of  alternating  vertical  and  horizontal  polarizations.  The  spatio-polarimetric 
correlation  coefficients  and  the  DOA  estimation  are  used  to  assess  the  performance  of  the 
array  reconfigurations. 

II.  Signal  Model 

Assume  L  source  signals  si(t),l  =  1, ...,  L,  are  incident  on  an  array  of  N  dual-polarized 
antenna  sensors.  We  use  M  and  ^  to  represent  the  orthogonal  vertical  and  horizontal 
polarizations  at  each  sensor.  The  received  data  for  each  polarization  is  the  linear  combi¬ 
nation  of  the  same  source  polarization  components  and  noise.  That  is,  the  signal  received 
at  the  nth  sensor  of  polarization  i,  i  =  v  or  h,  is 

Xn  (t)  =  53  anA] (t)  +  Un  (*) ,  (!) 

1=1 

where  a^(  represents  the  propagation  coefficient  of  the  Ith  source  with  polarization  i,  (  t) , 
to  the  nth  array  sensor,  and  r${t)  is  the  noise  component  at  the  same  polarization. 

The  data  vector  at  each  polarization,  xW(t),  i  =  v  or  h,  is  expressed  as 

xM(t)  =yW(i)+nM(t) 

=  AWSM(i)+nM(*)  (2) 

=  AWpMs(*)  +  nW(t), 

where  y W  (t)  and  nW  (t)  are,  respectively,  the  noise-free  data  vector  and  the  noise  vector, 
AW  =  [aP,...,ag]  is  the  data  mixing  matrix,  and  sW(t)  is  a  vector  representing  the 
ith  polarization  components  of  the  source  signal  vector,  s (t).  In  addition,  the  polarization 
information  can  be  decoupled  from  the  source  waveform  and  separately  captured  in  matrix 
pH  =  diag \pi\  ■  ■  ■  ,p^],  where 

pf]  —  cos  (71)  and  (3) 

are  the  coefficients  of  the  Zth  source  projected  on  the  v  and  h  polarizations.  In  (3), 

7 1  G  [0,  |]  is  the  polarization  angle  that  determines  the  magnitude  ratio  between  the 
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two  polarizations,  whereas  r]i  €  (—7 r,  i r]  is  the  phase  difference  between  the  two  polariza¬ 
tions.  In  this  paper  we  only  consider  linearly  polarized  signals,  i.e.,  rft  =  0,  l  =  1, . . . ,  L. 
Without  loss  of  generality,  it  is  assumed  that  the  norm  of  the  spatial  signature  vector 
41  =  [a\\,  •  •  •  ,  a$z]r  for  each  source’s  polarization  is  N,  where  l  =  1,  •  •  •  ,  L,  i  =  v,  h , 
and  T  denotes  transpose.  From  (3),  it  is  clear  that  p i  =  [pf1  pf^]T  has  a  unit  norm.  Thus, 
the  strength  of  the  source  signals  is  absorbed  in  the  magnitude  of  si(t). 

III.  Comparison  between  different  array  configurations 

We  now  compare  an  iV-sensor  dual-polarization  linear  array  with  two  configurations  of 
2JV-sensor  dual-polarized  arrays,  each  is  composed  of  single-polarized  elements,  as  shown 
in  Fig  V.  Each  antenna  of  the  dual-polarized  array  of  dual-feed  sensors  has  dual-feed 
points,  as  shown  in  Fig.  1(c).  In  the  first  dual-polarized  array  of  single-polarized  elements 
configuration  (Configuration  A  in  Fig.  1(a)),  the  array  consists  of  N  vertically  polarized 
sensors  followed  by  N  horizontally  polarized  sensors,  whereas  in  the  second  dual-polarized 
array  of  single-polarized  elements  configuration  (Configuration  B  in  Fig.  1(b)),  the  array 
consists  of  2 N  sensors  alternating  between  horizontal  and  vertical  polarizations. 


The  combined  data  vector  received  at  a  dual-polarized  antenna  array  is  expressed  in  the 
following  vector  format 


x(t) 


xM(£) 

AMpM 

s(t)  + 

n  M(t) 

xM(t) 

A^ipW 

nW(t) 

(4) 


The  noise  elements  are  modelled  as  stationary  and  white  complex  Gaussian  processes  with 
zero  mean  and  variance  a2  in  each  polarization,  i.e., 


E  [n(i  +  r)nff(t)]  =  a2S(r) l2N, 


(5) 


where  5(r)  is  the  Kronecker  delta  and  I 2n  denotes  the  2 N  x  2N  identity  matrix. 

In  the  underlying  dual-polarized  antenna  arrays,  we  consider  the  closeness  and  distinc¬ 
tions  of  two  sources  using  the  correlation  in  both  the  spatial  and  polarimetric  domains.  It 
is  clear  from  (4)  that 
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represents  the  joint  spatio-polarimetric  signature  of  the  Zth  signal,  l  —  1,  •  •  •  ,  L. 

A.  Dual- Polarized  Array 

We  consider  a  structured  mixing  matrix  for  convenience.  In  this  case,  aj^  =  a}^  =  a;. 
Accordingly,  the  joint  spatio-polarimetric  signature  of  the  Zth  signal  for  a  dual-polarized 
array  can  be  expressed  as 

a \D)  =  a i®  pjD),  (7) 

where  ®  denotes  the  Kronecker  product  operator.  In  conventional  single-polarized  antenna 
arrays,  a  high  spatial  correlation  may  prohibit  source  resolution  [7].  Therefore,  reduction 
of  the  source  spatial  correlations  becomes  important  which  can  be  achieved  by  changing 
the  array  manifold.  The  new  manifold  is  comprised  of  both  source  spatial  and  polarization 
signatures.  The  spatio-polarimetric  correlation  coefficient  is  defined  as  [8] 


where  H  represents  the  Hermitian  and  p\°2  =  [pm  ^  p[D)  is  the  polarimetric  correlation 
coefficient  between  sources  l  and  m,  and  (dl>m  —  (l/Aja^a*  is  the  spatial  correlation 
between  these  sources  for  an  ZV-element  linear  array.  The  iV-element  array  response  vector, 
a /,  is  expressed  as 

a/  =  [1  e~j9‘  ■  ■  ■  e-W-W'f,  (9) 

where  Oi  =  2ir(d/\)  sin  </)(,  d  is  the  interelement  spacing  of  the  array,  A  is  the  wavelength 
and  (j>i  is  the  DOA  of  the  Zth  source.  The  polarimetric  correlation  coefficient  between  two 
sources  is  given  by 

Pi,m  =  cos  H  cos  Tm  +  sin  7 i  sin  jm  =  cos (7,  =p  7m)  •  (10) 

B.  Configuration  A 

When  using  a  dual-polarized  array  composed  of  the  augmentation  of  N  vertically  po¬ 
larized  sensors  and  N  horizontally  polarized  sensors,  the  array  response  is 

a.\A)  =a,®p{A),  (11) 
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with 


(A) 

Pi  = 


cos  7 1  sin  7 ie 


-jN6, 


(12) 


Resembling  (8),  the  spatio-polarimetric  coefficient  is 

(13) 

where  p\A^  is  a  function  of  the  polarimetric  correlation  and  DOAs  and  is  expressed  as 

C  =  (pLA))%!A) 


-jN(e,-em) 


cos  7 1  cos  7m  +  sin  ji  sin  jme 


The  spatio-polarimetric  correlation  is  reduced  for  p\A^  <  1. 

C.  Configuration  B 

Similarly  the  spatial-polarimetric  array  response  for  configuration  B  is 

a[S)  =  (a*  0  a*)  <g>  pjB), 


(14) 


(15) 


where  ©  denotes  the  Hadamard  product.  The  element-by-element  product  of  the  N- 
element  array  response,  a;  results  in  the  array  response  vector  of  a  decimated  array  with 
twice  the  interelement  spacing  as  the  original,  and 

iT 


-j6i 


(16) 


( b )  r 

Pi  =  |c0S7i  sm 7ie~ 

Subsequently  the  spatio-polarimetric  coefficient  is 

C = (i7) 

where  pt  'm  =  ^  (a^  ©  am)  (a;  ©  a()  is  the  spatial  correlation  for  a  A-element  decimated 
array.  The  polarimetric  correlation,  which  also  contains  a  spatial  component,  is 

P{iB2  =  (Pm])HP{iB) 

=  cos  7 /  cos  7 m  +  sin  sin  7 (18) 

It  is  evident  from  (17)  that  the  array  used  in  configuration  B  with  half  a  wavelength 
interelement  spacing  generates  a  grating  lobe  problem  when  the  spatial  separation  between 
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the  two  sources  is  high.  To  avoid  the  grating  lobes,  the  array  in  configuration  B  must  be 
spatially  over-sampled  such  that  d  <  j.  If  the  interelement  spacing  is  reduced  to  d  = 
then  the  aperture  of  the  array  in  configuration  B  becomes  —4-  which  is  roughly  the  same 
as  — -  for  the  dual-polarized  array. 

IV.  Simulations 

We  compare  the  performance  of  three  array  configurations  for  polarimetric  MUSIC- 
based  DOA  estimation.  We  also  relate  their  performance  to  the  sources’  spatio-polarimetric 
correlation.  Figure  2  compares  the  spatial  correlation  of  an  N  =  4  element  array  with 
d  =  A/2  and  a  decimated  array  with  d  =  A  for  a  source  with  a  fa  =  10°.  The  effect  of  the 
spatial  under-sampling  is  clearly  demonstrated. 

Figure  3  compares  p{^,  p[^  and  p[B-]  for  (pi  —  10°  and  7i  =  45°.  It  can  be  seen  that 
p[°2  is  independent  of  72  meanwhile  p'^}  and  p\B]  vary  with  changes  in  72.  A  grating 
area  is  shown  in  Figs.  3(b)-(d)  due  to  the  exponent  in  (14)  and  (18).  Additionally,  spatial 
over-sampling  for  the  array  in  configuration  B  reduces  the  grating  area  in  Fig.  3(d)  as 
compared  to  Fig.  3(c).  It  is  worth  noting  that  the  same  spatial  over-sampling  eliminates 
the  grating  problem  in  the  spatial  correlation, 

Figure  4  depicts  the  spatio-polarimetric  correlations  and  /3^  versus  the  DOA 

and  polarization  angle  parameters,  fa  and  72.  The  dual  polarized  array  does  not  exhibit 
a  grating  problem  as  can  be  seen  in  Fig.  4(a).  The  array  in  configuration  A  shows  some 
grating  area  due  to  the  exponent  in  (16),  however,  the  grating  area  is  lower  than  that  in 
Fig.  3(b)  due  to  the  inclusion  of  the  spatial  correlation.  The  effect  of  spatial  over-sampling 
is  evident  in  the  reduction  of  the  grating  area  in  Fig.  4(d),  as  compared  to  that  of  the 
decimated  array  exhibited  in  Fig.  4(c). 

Figure  5  shows  the  performance  comparison  between  the  three  arrays  at  hand  for  two 
different  apertures.  In  Fig.  5(a),  the  apertures  of  all  the  arrays  is  kept  equal  to  the  aperture 
of  the  dual-polarized  array.  The  performance  of  the  dual-polarized  array  is  comparable  to 
that  of  the  array  in  configuration  B.  In  this  case,  the  array  in  configuration  B  does  exhibit 
the  grating  lobe  problem,  meanwhile  the  array  in  configuration  A  is  spatially  over  sampled. 
In  Fig.  5(b),  the  apertures  of  all  the  arrays  is  the  same  as  the  aperture  of  the  array  in 
configuration  A.  The  performance  of  the  dual-polarized  array  is  also  comparable  to  the 
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array  in  configuration  B.  On  the  other  hand,  the  performance  of  the  array  in  configuration 
A  suffers  from  further  loss  of  resolution  compared  to  the  other  two  arrays  since  its  spatio- 
polarimetric  correlation  is  not  as  low  as  its  counterparts.  It  is  worth  noting  that  the 
dual-polarized  array  and  the  array  of  configuration  B  both  suffer  from  the  grating  lobe 
problem  in  Fig.  5(b). 

V.  Conclusion 

The  paper  compared  the  DOA  performance  of  different  configurations  of  multi-antenna 
receivers  having  single-polarized  and  dual-polarized  double-feed  antennas.  All  arrays  have 
the  same  number  of  antennas  with  one  half  of  the  antennas  of  vertical  polarizations  and 
the  other  half  of  horizontal  polarizations.  The  DOA  performance  of  each  array  is  viewed  in 
terms  of  the  respective  polarimetric  spatial  correlation.  To  establish  a  common  reference, 
the  correlation  is  expressed  in  terms  of  a  product  of  spatial  correlation  of  a  reference  array 
and  a  modification  that  reflects  the  change  in  array  configuration.  The  paper  contrasted 
the  reduction  in  correlation  with  the  grating  lobes  that  arise  due  to  using  the  same  array 
span. 
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(a)Configuration  A 


(b)  Configuration  B 


(c)Dual-polarized  array 
Fig.  1.  Array  configurations 
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MICRO-DOPPLER  EXPERIMENTS 

Pawan  Setlur,  Moeness  Amin,  and  Fauzia  Ahmad 

Abstract 

Targets  in  real  life  have  vibrating  and  rotating  parts  in  addition  to  the  translations  endured  by 
the  target.  The  target  itself  might  either  be  rotating  or  vibrating  with  no  translation,  a  motion 
is  described  mechanically  as  simple-harmonic  motion  (SHM).  Such  nonlinear  motions  then 
induce  additional  frequency  modulations  around  the  Doppler,  this  is  termed  micro-Doppler 
effect.  Micro-Doppler  has  many  applications  for  target  recognition  and  classification.  We 
have  conducted  several  indoor  experiments  in  the  Radar  Imaging  Lab  at  the  Center  for 
Advanced  Communications  at  Villanova  University,  for  detecting  micro-Doppler. 
Experiments  were  designed  to  image  targets  which  induce  micro-Doppler  like  signatures. 
The  targets  were  primarily  a  rotating  fan  and  a  toy  train  moving  on  a  circular  track.  The 
experiments  were  conducted  using  CW  (Continuous  Wave)  radars.  These  radars  measure 
only  the  real  component  of  target  returns.  Time-frequency  techniques  were  used  to  depict  the 
instantaneous  frequency  of  the  target.  The  objective  of  the  experiments  was  to  compute  the 
frequency  of  vibrations  or  rotations.  In  the  ensuing  sections,  we  introduce  the  problem 
statement,  and  propose  a  model  for  Micro-Doppler.  Experimental  setup  is  described,  and 
simulation  and  experimental  results  are  also  presented. 
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I.  INTRODUCTION 


The  effect  of  nonuniform,  rotational  or  non-linear  motion  on  target  signatures  and 
characteristics  was  extensively  studied  in  the  past  [1],  However,  only  recently,  returns  from 
rotational/vibrational  motions  of  targets,  termed  micro-Doppler,  have  been  investigated  [2,  3], 
The  term  “micro”  is  primarily  used  for  two  main  reasons  i)  Smaller  sidebands  about  the 
Doppler,  ii)  Micro  motions  of  targets,  compared  to  large  displacements  involved  in 
translations  [4].  Micro-Doppler  captures  certain  dynamics  of  the  target,  which  aid  in  target 
classification  and  recognition  [3,  4,  5],  as  every  target  has  distinctive  micro  motion  dynamics. 
Fourier  based  methods  have  inherent  resolution  problems.  Hence,  high  resolution  time- 
frequency  methods  are  employed  [3,  4]. 

Vibrational/rotational  target  endure  SHM.  Hence,  irrespective  of  the  model  chosen, 
whether  a  rotation  model  or  a  vibrational  model,  the  target  signature  will  be  sinusoidal  in  the 
time  frequency  plane.  Moreover,  the  Doppler  frequency  also  exhibits  SHM,  with  the  same 
resonant  frequency  as  the  target. 

II.  MODEL 

A  micro-Doppler  signal  can  be  parameterized  by, 

x(t)  =  exp(y'2^cos(®v/  -  ^))  (1) 

where  cov  is  the  frequency  of  vibration/rotation. 

The  instantaneous  frequency  IF(l)  can  be  written  as, 

IF(t)  =  -ov  sin(©v/  -  (f)  (2) 

which  satisfies 

£ r(lF(t))+a;iF(t)=0  (3) 

at 
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Equation  (3)  is  the  differential  equation  for  undamped  SHM  (Simple  Harmonic  Motion). 
Hence,  the  instantaneous  frequency  IF(t)  exhibits  SHM  with  the  same  resonant  frequency.  It 
must  be  noted  that  rotation  is  also  SHM.  This  can  be  proven  when  we  consider  the  motion  to 
be  projected  onto  the  diameter  of  the  circular  trajectory  of  the  target.  Therefore,  in 
practicality,  we  do  not  need  another  rotation  model  to  describe  targets  having  turning 
motions.  Figure  1  is  the  3D  schematic  of  a  radar  and  a  target  enduring  SHM.  If  the  target  has 


an  instantaneous  range  given  by 

r(t)  =  dcos(<yv/)  (4) 

d  =  d[cos(ttr  )cos(/?, ),  cos(er,  )sin(/?, ),  sin(«,  )F  (5) 

Then,  the  instantaneous  frequency  of  such  a  setup  can  be  written  as 

I?'/  avd  sm(o)vtXnr)T  nt  (6) 

nr  =  (cos(ar  )cos(flr ),  cos(ar  )sin(/?r  ),sin(a,))r  (7) 

nt  =  (cos(a,  )cos(/?, ),  cos(«r,  )sin(/I, ),  sin  (or,  F  (8) 


Here,  eov  is  the  resonant  or  fundamental  frequency  of  the  SHM  of  the  target,  d  is  the 

maximum  displacement  of  the  target  from  the  origin,  and  fc  is  the  carrier  frequency  of  the 

radar.  The  pairs  ( a, ,  (5t )  and  (ar ,  pr )  are  the  elevation  and  azimuth  of  the  target  and  radar 

respectively.  Clearly,  equation  (6)  gives  rise  to  a  sinusoidal  signature  in  the  time-frequency 
plane. 
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Figure  1.  Geometry  of  the  scene 

III.  HARDWARE  SETUP 

The  24  GHz  and  10  GHz  CW  (Continuous  wave)  radars  were  used  for  imaging  targets 
undergoing  micro-Doppler  like  motion.  The  radar  units  used  measure  only  the  real 
components  of  the  baseband  returns.  Hence,  additional  signal  processing  must  be  performed 
to  properly  analyze  the  returns.  An  A/D  converter  was  used  to  digitize  the  analog  returns. 
The  A/D  converter  can  operate  on  8  different  single  channels  simultaneously.  An  application, 
written  in  Visual  Basic®  on  Microsoft  .Net®  platform,  is  used  to  set  the  sampling  rate  and 
communicate  with  the  A/D  converter.  Figure  2  depicts  the  collection  of  main  hardware  and 
radars.  Figure  3  depicts  a  snapshot  of  the  application  used  to  interface  the  A/D  with  a  laptop. 
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Figure  3.  Snapshot  of  application. 


We  chose  CW  radars  as  they  are  affordable,  robust  to  clutter,  and  perform  exceptionally  well 
in  highly  cluttered  environments  frequently  encountered  in  indoor  imaging. 

IV.  DATA  COLLECTION 

We  have  performed  three  experiments  producing  various  types  of  micro-Doppler  signatures. 
These  are:  high  micro-Doppler  for  rotating  fans,  very  low  micro-Doppler  for  a  toy  train 
moving  on  a  circular  track,  and  low  frequency  components  for  a  comer  reflector  moving 
back  and  forth  at  a  slow  rate.  We  will  describe  the  experimental  setup  in  the  ensuing  sub¬ 
sections.  However,  we  digress  to  bring  other  important  aspects  to  light.  If  x(t)  is  the  real 
signal,  and  H{}  denotes  the  Hilbert  transform  operation,  the  analytic  signal  y(t)  is  then 
given  by 

y(t)  =  x(t)+jH{x(t)}  (9) 

The  baseband  micro-Doppler  signal  has  no  carrier  and  has  an  instantaneous  frequency  on 
both  sides  of  0  Hz.  Hence,  the  transformation  in  (9)  will  mirror  the  negative  instantaneous 
frequency  onto  the  positive  side.  This  can  be  seen  in  Figure  4,  which  shows  the  pseudo 
Wigner-Ville  distribution  of  the  analytic  signal,  after  performing  transformation  (9)  on  a 
simulated  real  micro-Doppler  signal.  Figure  5  is  the  pseudo  Wigner-Ville  distribution  of  the 
complex  signal.  Note  that  Figure  5  is  the  actual  target  signature. 
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Figure  4.  Analytic  signal’s  signature. 


Figure  5.  Complex  signal’s  signature  (Actual  Target  signature). 
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The  inherent  drawback  of  using  real  signals,  therefore,  is  that  the  frequency  components  are 
symmetric  (replicated)  on  either  side  of  the  frequency  axis.  In  our  case,  since  the  signal 
frequency  spectrum  is  centered  at  0Hz,  we  expect  to  see  replicas  around  0Hz.  For  micro- 
Doppler  detection,  this  replica  can  be  misinterpreted  as  another  target.  For  example,  in  Figure 
6,  we  compute  the  pseudo  Wigner-Ville  for  a  real  micro-Doppler  signal  with  a>v  =  In  Hz 
and  (f>  =0  in  (1).  The  Figure  indicates  two  targets  having  vibrations  180°  out  of  phase  with 
each  other.  In  reality,  however,  there  is  only  one  target.  Since  we  have  apriori  information 
about  the  experiments,  radar  output,  and  target  positions,  we  can  ignore  the  replica  and 
consider  only  the  target  signature. 


Figure  6.  Real  signal’s  signature. 
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V.  EXPERIMENTAL  SETUP 

We  have  conducted  three  different  experiments,  showing  the  micro-Doppler  signatures  of 
objects  that  could  be  found  in  indoor  radar  imaging,  namely  (a)  Comer  reflector  moving  back 
and  forth  over  an  arms  length,  (b)  Rotating  fan,  (c)  Toy  train. 

V.l.  Corner  Reflector 

The  comer  reflector  echoes  most  of  the  transmitted  energy  back  to  the  radar.  The  comer 
reflector,  used  in  the  experiment  and  shown  in  Figure  7,  was  held  by  a  person  and  moved 
back  and  forth  slowly  in  front  of  the  radar.  The  initial  distance  from  the  radar  is  5  ft.  The 
sampling  frequency  used  is  100Hz  for  a  duration  of  5  seconds.  The  10  GHz  radar  was  used  to 
image  the  target. 


Figure  7.  Corner  reflector  used  in  experiment  A. 
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V.2.  Rotating  fan 

The  fan  is  portable  and  has  a  height  of  approximately  3  feet  from  the  ground.  It  has  3  plastic 
blades.  We  covered  one  blade  with  aluminum  foil  so  that  the  returns  from  the  other  two 
blades  are  not  seen  by  the  radar.  The  fan  and  the  experimental  setup  are  shown  in  Figure  8. 
High  frequency  micro-Doppler  components  are  expected  from  a  rotating  fan.  The  sampling 
frequency  must  be  high  enough,  so  that  the  window  used  in  the  time-  frequency  distribution 
can  capture  the  non-stationarity  of  the  signal.  We  simulated  a  rotating  fan  using  <yv  =  10Hz 
and  ^  =  0  in  equation  (1).  We  set  the  sampling  frequency  to  100Hz  and  1000Hz,  and 
computed  the  pseudo  Wigner-Ville  distribution  for  both  cases.  In  Figure  9,  the  fan’s  returns 
are  sampled  at  100Hz  and  the  signal  length  is  1000  samples.  In  Figure  10,  the  sampling 
frequency  is  1000Hz  and  number  of  samples  are  1000.  The  window  size  for  the  pseudo 
Wigner-Ville  is  the  same  in  both  cases.  Figure  9  does  not  depict  the  signal’s  instantaneous 
characteristics,  but  rather  the  overall  spectrum  of  the  harmonics.  On  the  other  hand.  Figure 
10  clearly  shows  the  instantaneous  frequency. 

The  sampling  frequency  used  in  the  fan  experiment  was  thereby  set  to  1000Hz  for  a 
duration  of  1  second.  The  base  of  the  fan  is  at  a  distance  of  4  feet  from  the  radar.  The  24  GHz 
radar  was  used  in  the  experiment. 

V.3.  Toy  train  on  a  circular  track 

The  toy  train  in  the  experimental  setup  can  be  seen  in  Fig.  11.  The  diameter  of  the  circular 
track  is  15.5  inches.  The  body  of  the  train  is  made  of  plastic.  The  returns  we  expect  are  low 
frequency  micro-Doppler  signals.  The  sampling  frequency  used  in  this  experiment  is  100Hz 
for  a  duration  of  1 7  seconds.  Aluminum  foil  was  draped  onto  the  body  of  the  train,  except  the 
wheels.  The  10  GHz  radar  was  used  for  this  experiment. 
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Normalized  frequency,  Hz 


Figure  8.  Layout  of  the  fan  experiment.  The  24  GHz  radar  can  be  seen  at  the  left  end  of  the 

table. 

Simulated  fan.  sampling  frequency,  100  Hz 
5 

4 

3 
2 
1 
0 
-1 
-2 
-3 
-4 

Figure  9.  Simulated  fan-sampling  frequency  of  100  Hz 


Normalized  frequency,  Hz 


Figure  11.  Layout  of  the  toy  train  experiment.  The  10  GHz  radar  can  be  seen  on  top  of  the  24 

GHz  radar. 
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RESULTS 


In  all  experiments,  we  used  the  pseudo  Wigner-Ville  distribution  to  illustrate  the 
instantaneous  frequency.  Figure  12  depicts  the  sinusoidal  signature  for  Experiment  (A).  Due 
to  the  replicas  and  the  aliasing  effects  in  the  time  frequency  plane,  the  time  period  cannot  be 
definitively  estimated.  The  results  of  Experiment  (B)  are  shown  in  Fig.  13,  where  we  can 
observe  two  signatures.  This  is  due  to  using  the  real  signal,  as  previously  discussed.  However, 
since  we  are  interested  in  only  estimating  the  resonant  frequency,  we  can  ignore  one  of  the 
signatures.  The  rotation  rate  of  the  fan  is  estimated  to  be  approximately  750rpm.  An 
interesting  observation,  from  0.6-1  second  in  Figure  13,  is  that  the  time  period  of  rotations 
decreases.  This  occurred  because  the  aluminum  foil  on  the  blade  unexpectedly  slackened  and 
hindered  the  smooth  operation  of  the  fan. 

Figure  14  illustrates  the  results  of  Experiment  (C).  As  observed  from  the  figure,  the 
toy  train  has  a  time  period  of  about  7.5  seconds  for  a  single  rotation.  A  stop  watch  was  also 
used  to  measure  the  time  period.  The  experimental  values  were  in  close  agreement  with  the 
measured  values.  At  time  instants  3.5  and  11  seconds,  pulse  like  events  (more  specifically, 
narrow  sine  like  functions)  are  seen  in  Fig.  14.  The  difference  in  time  for  the  events  is  equal 
to  the  time  period  of  one  full  rotation.  The  train  has  a  metallic  rod  connecting  the  rear  wheels, 
one  end  of  which  is  highly  reflective.  The  pulse  like  events  are  caused  by  the  highly 
reflective  end  of  the  rod,  which  is  exposed  to  the  radar  over  a  short  period  within  one  rotation. 
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Figure  13.  Results  of  experiment  B 
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Figure  14.  Results  of  experiment  C 
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